An illustrative example of the relationship between dissipation and relative entropy 
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PACS numbers: 05.70.Ln, 05.40.-a 



I. INTRODUCTION 

During a thermodynamic process in which a system, in 
contact with a thermal reservoir, evolves from one state of 
thermal equilibrium A, to another B, the average work 
performed on the system must exceed the free energy 
difference: (W) > AF = F B - F A 1\. The work in 
excess of AF, i.e. the dissipated work (W<ii SS ) = (W) — 
A.F, quantifies the amount of energy irretrievably lost to 
the surrounding thermal environment. Recently, Kawai, 
Parrondo, and Van den Broeck Q have related (Wdiss) to 
another measure of irreversibility, roughly speaking the 
distinction between the forward and reverse directions of 
the "arrow of time" . Specifically, they have shown that 



f3(W diss ) >D( PF \\p R ), 



(1) 



where pf and p^ are time-dependent phase-space den- 
sities describing the evolution of the system from A 
to B, and from B to A, respectively; and D{p\\p') — 
J p\n(p/p') denotes the relative entropy 0, a measure of 
the distinguishability between two distributions. When 
the system in question evolves detcrministically, under 
Hamilton's equations, then Eq.[T]is an equality. However, 
if the system description is coarse-grained, or explicitly 
stochastic, then the relative entropy only provides a lower 
bound on the dissipated work. 

Connections between dissipation and temporal asym- 
metry similar to Eq.Q]have also been established by other 
authors. Maes [4[ and later Maes and Netocny Q have 
obtained a correspondence between thermodynamic mea- 
sures of dissipation and relative entropies between for- 
ward and reverse distributions in path space; and Gas- 
pard [|| has connected this result to the dynamical ran- 
domness that characterizes nonequilibrium steady states. 
In the context described in the previous paragraph, one 
of us [7| has derived a relation analogous to Eq. [I] but 
expressed in terms of distributions in path space (see 
Eq.EJ). 

The goal of the present paper is to illustrate Eq.Q]using 
the simple model of a particle inside a one-dimensional, 
moving harmonic well. This model is both analytically 
tractable and experimentally relevant [1, H, EH E2, EE 
E3| ■ After briefly reviewing the central result of Ref. [2| 



in Section |TT1 we introduce our model at the beginning 
of Section IIIIl We then solve the model and illustrate 
Eq. [1] when the system evolves under Hamilton's equa- 
tions (Section IIII A[) . and both overdamped and iner- 
tial Langevin dynamics (Sections IIII Bl and IIII CI respec- 
tively). In Section llVI we extend our analysis, and show 
that the bound in Eq. [1] can be improved by specifying 
the phase space density at two times, rather than one. 
Gomez-Marin et al. [15|, [l6| have also studied this model 
in the context of Eq. [1] and in Section IIVI we briefly dis- 
cuss the relationship between our results and theirs. 



II. BACKGROUND 

Consider a classical system with N degrees of freedom, 
described by the coordinates x = {xi, . . . , xn} and con- 
jugate momenta p = {pi, . . . ,Pn}- Let z — (x, p) de- 
note a point in phase space, and let H(z, A) denote a 
parameter-dependent Hamiltonian for this system. We 
assume that this Hamiltonian is time-reversal invariant: 
H(z*, A) = H(z, A), where the asterisk denotes the rever- 
sal of momenta, p — * —p. 

The term thermodynamic process will indicate a se- 
quence of events whereby the system evolves in phase 
space as the external parameter A, is varied according 
to an arbitrary schedule, or protocol X t (also labeled as 
X(t)), from an initial time t = to a final time t = r. 
During this interval of time, the microscopic evolution of 
the system is specified by a trajectory z tl or z(t), and 
p{z, t) will denote the time-dependent phase space den- 
sity describing an ensemble of such trajectories. As in 
Ref. [|J we will explicitly consider two such processes: 
one defined by a forward protocol Xf, during which the 
parameter is varied from an initial value A^ = A to a fi- 
nal value A^T = B, the other defined by a reverse protocol 
X* from Xn — B to X?- = A. 



F t , from Xq = B to A« 
Prior to the start of either the forward or the reverse 
process (t < 0), the system is brought to equilibrium by 
weak contact with a thermal reservoir at inverse temper- 
ature (3. As a result, the initial phase space density is a 
canonical distribution at A = A or B. The system might 
subsequently be thermally isolated from the reservoir, or 
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else it might remain in contact with the reservoir. In the 
former case, Hamilton's equations govern its evolution 
from t = Q to t = r, while in the latter case we will use 
stochastic dynamics to model the random effects of the 
environment. In either situation the work performed on 
the system during this interval of time is given by the 
following functional of the trajectory: 



w[zt] = dtx—(z t ,\ t ). 



(2) 



For the case of Hamiltonian dynamics, the average dis- 
sipated work can be expressed as (see Ref. @ for the 
details of the derivation) 



P(W diss ) = J dzp F (z,t) In 
= D (p F (z,t)\\p R (z 



PF(z,t) 

Pr(z*,t- t) 
\r-t)). 



(3) 



Since relative entropy is a measure of the distinguishabil- 
ity of two probability distributions [|[ , Eq. [3] relates the 
dissipation of energy to the ease with which a process can 
be distinguished from its time-reversal @, [H, [H, [l7| • 

Eq. [3J is an equality because the underlying Hamilto- 
nian dynamics are deterministic. If the system remains in 
contact with the reservoir during the process, we instead 
obtain an inequality, Eq. [TJ There are two ways to argue 
this, both of which make use of a key property of the 
relative entropy between two distributions, namely that 
it decreases when the distributions are projected onto a 
smaller set of variables [3[. When the system remains in 
contact with a reservoir, then we can view the system 
and the reservoir as two sub-systems that together form 
a large, isolated "super-system" to which Eq. [3J applies: 
the dissipated work gives the relative entropy between 
the two distributions in the full phase space. Upon pro- 
jecting out the reservoir variables, the relative entropy 
decreases, and we obtain Eq. [TJ as discussed in Ref. Q. 
Alternatively, we note that an equality analogous to Eq.[3] 
can be formulated for distributions in path space, rather 
than phase space Q: 

(3(W diss ) = DiVF^WR^)). (4) 

Here the trajectory j F describes the evolution of the sys- 
tem during a given realization of the forward process, and 
Pf{i F ) is the probability distribution of such trajecto- 
ries; 7^ and Vr are defined in a similar manner for the 
reverse process. Equation 2] holds for both deterministic 
and stochastic evolution. If we now project from path 
space onto phase space, e.g. V F (j F ) — > p F (z 7 t), then 
the relative entropy decreases, and we again obtain the 
inequality, Eq[TJ (References [HI, 03 refer to this projec- 
tion procedure as "coarse-graining in time".) In either 
case, the decrease of relative entropy has a simple in- 
terpretation: when we discard microscopic information - 
either about the state of the reservoir, or about states 
of the system itself at times other than the specified in- 
stant, t then we diminish our ability to distinguish be- 
tween the forward and the reverse process. This sug- 
gests that the more microscopic information we retain, 



the closer the value of D is to its upper bound /3(Wdi ss )- 
We investigate this in Section HVi where we consider a 
generalization of Eq. [1] in which the microstate of the 
system is specified at two instants in time, rather than 
one. References @, [TH, [IH, [TtJ contain a related analysis 
specifically addressing the loss of microscopic informa- 
tion during coarse-graining and how this loss affects the 
value of the relative entropy. 



III. PARTICLE IN A MOVING HARMONIC 
WELL 

In this section we illustrate Eq. [1] through an explicit 
calculation of the average dissipated work and relative 
entropy for a particle in a one-dimensional, moving har- 
monic well. We begin by specifying the model and ob- 
taining useful preliminary results. Then in section [III Al 
we evaluate the average dissipated work and relative en- 
tropy for a system following Hamiltonian dynamics. In 
sections MI Bl and MI CI we repeat the calculation with 
the system modeled by Langevin dynamics, in different 
limiting regimes. 

Our system is a particle of mass m trapped in a har- 
monic well with spring constant k: 



H(x,p, A) 



P 

2m 



+ --{x— A) 2 



(5) 



We will consider processes during which the center of 
the well is moved cither rightward or leftward at con- 
stant speed u. These correspond to forward and reverse 
protocols, X F (t) — ut and Xa(t) — u(t — t). Explicit 
expressions for the initial equilibrium densities are given 
by the following Gaussians: 



Pf(z,0) 



Pr(z,0) 



JL 

A 
2tt 



m 
~k 



exp 



exp 



-0 



kx 2 



V 

2m 

P 2 , k 2 

[X — UT) 

2m 2 y ' 



(6a) 



(6b) 



For this system AF — by translational symmetry, 
therefore for the forward process 



(Wdiss) = -uk / dt[x F (t) - ut], 
Jo 



(7) 



where x F (t) is the average position during the forward 
process. 

Since the time-dependent densities p F {z,t) and 
Pr(z, t) will prove to be Gaussians for all the cases con- 
sidered in this paper (see Sections MI Al - MI C[) , it is 
useful here to establish uniform notation for the descrip- 
tion of two-dimensional Gaussian distributions fa( z ) — 
fc{x,p). Such a distribution is uniquely determined by 
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the moments (means and covariances) 



X = 


p 

j dzf G (z)x, 




(8a) 


V = 


r 

/ dzf G (z)p, 




(8b) 


°l = 


j dzf G (z)(x 




(8c) 


°l = _ 


j dzf G (z)(p 


-pf, 


(8d) 


(7 X p — 


j dzf G (z)(x 


- x)(p-p). 


(8e) 



An explicit expression for f G (z) in terms of these mo- 
ments is 



fa(z) 



1 



27r\/det a 



exp 



~(z -zf -a 1 -{z- z) 



, (9) 



where z is a vector in phase space and a is the covariance 
matrix: 



a x a xp 



(10) 



Lastly, the relative entropy between two Gaussian distri- 
butions, /g(z) and ga(z*), is (see Eq. [3]) 



D(f G (z)\\g G (z*))=-l + 



hi 



/ det tj g 
\det 07 



+ o( z / ~ z g) T ' a g 1 ' ( z / - z s)> 



(11) 



where er*„ = ~u xp and all other elements of a* are unal- 
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and the initial variances 



«0) = 0, (13b) 



which are the same for the forward and reverse processes. 
Solving Eqs. [T^] leads to a set of linear equations for the 
positions and momenta in terms of their initial condi- 
tions. Combining these solutions with Eqs.[8land[T3lleads 
to the solutions 



x F (t) 



XR{t) 



sin(cji) 



sin(wt) 



PF(t) = mu [l — cos(cjt)l , 
(14a) 

PR.{t) = TOM[cOs(wi) — l] , 

(14b) 



= ^(*) = i» M«)=0, (14c) 



/3 



where w 2 = fc/m. 

We can now explicitly verify Eq. [1] The relative en- 
tropy between the forward phase space density pp(z,t) 
and the reverse phase space density pr(z*,t — t) is de- 
termined by plugging Eqs. Qj] into Eq. [Til 

D(p F (z,t)\\p R (z*,T- t)) = /W(l - cos(wt)). (15) 

Comparing this with 

(3(W diss ) = Pmu 2 (1 - cos(wr)) , (16) 



computed from Eqs. I14al and [71 we find the predicted 
result /3(Wdi ss ) = D(p F \\p R ). 



A. Hamiltonian Dynamics 

We now consider the case in which the system is ther- 
mally isolated from the reservoir after the initial equi- 
libration stage, and thereafter evolves under Hamilton's 
equations as the harmonic well is translated leftward or 
rightward, 

x F = pf/tti, pf = —k{xp — ut), (12a) 
xr = Pr/tu, pr — -k(xR - ut + ut). (12b) 

To solve for the evolution of the phase space densities 
Pf and pr, we observe that since the initial density is 
Gaussian (Eqs. O, and the equations of motion are lin- 
ear (E qs. [T2"|) . the distribution remains Gaussian for all 
times [19f . Thus, we need only determine the means and 
(co)variances as functions of time. From Eqs. [5] we have 
the initial means 

x F (0)=0, Pf(0) = 0, x r (0) = ut, p R (0)=0, 

(13a) 



B. Overdamped Langevin Dynamics 

We now imagine that the system remains in contact 
with the thermal reservoir throughout the process, and 
we will use Langevin dynamics to model the presence of 
the reservoir. As discussed in Section [Til Eq- UJ should 
apply as a strict inequality in this case. 

In this section we consider the overdamped limit, in 
which the momentum effectively equilibrates instanta- 
neously, and as a result the momentum does not con- 
tribute to the relative entropy. We therefore focus on the 
position dynamics. For the forward process the Fokker- 
Planck equation for pp{x,t) is 

d k d Id 2 

—p F (x,t) = -— [(x- ut)p F {x,t)} + ——p F ( x ,t), 

(17) 

where 7 is the friction coefficient. 

To solve Eq. [T7J we recognize that an initially Gaussian 
distribution will remain Gaussian for all time under the 
evolution of Eq. [TTJ as can be checked by substitution 
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(see for example [a [2y]) . Thus, Pf(%, i) is Gaussian with 
mean and variance, 



Xp it) = ut 



(l_ e -w/ 7 ) 



1 

~0k' 



(18) 



For the reverse process we replace t with r — t in Eq. 1171 
The solution is a Gaussian distribution with mean and 
variance, 



xr 



(t) = u(r-t) + ^(l-e-^-^), 



4 



1 

w 

(19) 

We now use these results to calculate the relative en- 
tropy and average dissipated work. Using Eq. [11] we ob- 
tain 



D(p F (x,t)\\p R (x,T-t)) 



2/3 7 2 u 2 



- fcT / 27 cosh 



k / T 



(20) 



The average dissipated work is obtained by substituting 
Eq. [THJinto Eq. [7] and evaluating the integral: 



P(W dh 



1 



(21) 



Verifying Eq. [T] requires demonstrating that the aver- 
age dissipated work is always greater than the relative en- 
tropy for any values of system parameters. To begin, we 
note that /3(Wdiss) does not depend on t and D(pf\\pr) 
obtains its maximum value D max at time t — t/2. Com- 
bining Eqs. 1201 and |2"T1 we get 



P(W diss ) - D{p F \\ PR ) > (3(W diss 
_ /3j 2 u 2 
~ k 
>0, 



C - (3 - e-^ 2 ) (l 



-C/2 



(22) 



where Q = kr/j is the scaled time. Using introductory 
calculus techniques, it is easy to verify that the bracketed 
quantity in Eq. [22] is non-negative for any ( > 0. As an 
example, we plot (3(W d i ss ) and D{p F \\p R ) in Fig. Q] as 
functions of time for fc, 7, u and [3 all set to one. 

In the quasi-static limit, namely u — *■ with ut fixed, 
both the dissipation and relative entropy approach zero, 
although at different rates. Specifically, from Eqs.l20"land 
|2T| /3(Wdi S s) ~ u and D(p F \\pn) ~ u 2 . This limiting be- 
havior can be justified on general grounds since we are 
considering continuous Markovian stochastic processes 
[2ll ] . A heuristic argument is as follows: consider such 
a continuous Markovian stochastic process perturbed by 
varying an external parameter A according to a specified 
protocol A(i). In the quasi-static limit, i.e. A ~ e <C 1, 
the phase space density is approximately 



(23) 



0.40 




FIG. 1: Comparison of ^(Wdiss) (dashed line) and D(pf\\pr) 
(solid line) in the overdamped limit. All parameters (fc, 7, u 
and f3) have been set to one in their respective units. 



where Pl s {z) is the unique stationary state with fixed 
external parameter A and 6P(z,t) is the first order cor- 
rection to the phase space density [22|]. Combining Eq. 
|2"31 with the definitions of average dissipated work (cf. 
Eq. [2]) and relative entropy, and taking the limit e — > 0, 
leads to /?(Wdi SS ) ~ e an d D{p F \\p R ) ~ e 2 . In the above 
model, A(f) = ut with e = u. 



C. Full Phase Space Langevin Dynamics 

In this section we analyze the same stochastic system 
without assuming the overdamped limit. (We do assume 
that the motion is not critically damped, i.e. 7 2 7^ Amk.) 
Using the same notation as the previous section pil Bp 
the Fokker-Planck equation for the phase space density 
of the forward process is 



d_ 

at 



p F {z,t) 



p d 

m dx 
d 



p F {z,t) 

k(x — ut) H Pf(z, t) 

m / 



dp 

-poV 2pF[Z ' 1 



(24) 



While Eqs. [24] can be solved exactly, the solution is 
complicated and unilluminating. An analytic solution is 
presented in Appendix [A] here, we illustrate Eq. [1] by 
plotting the ratio P{W d i SS ) / D(p F \\p R ) in Fig. Ofor var- 
ious values of the friction coefficient 7. Roughly speak- 
ing, 7 measures the coupling between the system and 
the reservoir: for larger 7, the reservoir and system in- 
teract more strongly; whereas, when 7 — > 0, the reservoir 
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Section IIII Bl The goal is to determine the pairwise 
probability P{x\,t\)Xo,to)- For the forward process, 
we can decompose the joint probability distribution 
into a conditional probability to be at x\ at t\ given 
the system was at xq at to, pf{x\, t\\XQ, to), and a 
one-time probability pF{xo,to), i-e. Pf{x\, t\\ Xo, to) — 
pF(xi,ti\xo,to)pp(xo,to). The expression for pp(xo,to) 
was determined in Section fill Bl (see Eq. IT5|) . The con- 
ditional probability is computed by solving the Fokker- 
Planck equation (see Eq. \T7\ with the initial condition 
8(x — xq) at to, i.e. the initial density is Gaussian with 
x(to) — xo and <J 2 {to) = 0. At time t\, the mean and 
variance for the forward process conditional probability 



x(t x ) =e- k ^-^^xo +u(h- toe- k{tl - to)h ^j 



11 (l _ e -fc(*i-*o)/7^ ) 



(25a) 



FIG. 2: Comparison of /3(Wdiss) / D(pf\\pr) for various values 
of the friction coefficient 7; the solid lines are, from bottom to 
top, 7 = 5 and 7 = 0.1. The dashed line represents Hamilto- 
nian dynamics (7 = 0) with /3(Wdiss) = D(pf\\pr)- All other 
system parameters (m, k, u and pT) have been set to one in 
their respective units. 



and the system decouple, with the result that the system 
evolves independently of the reservoir, under determinis- 
tic Hamiltonian dynamics. Figure [2] clearly shows that as 
7 decreases, the relative entropy approaches the average 
dissipated work. Decreasing 7 weakens the interaction 
between the system and the reservoir, resulting in less 
microscopic information leaking into the reservoir vari- 
ables. The result is more information is contained in the 
system variables, reflected by an increase in the relative 
entropy. The flow of microscopic information between 
system and reservoir degrees of freedom and its relation 
to dissipation has previously been discussed in the con- 
text of other models [H, \\% . 



IV. TWO-TIME PHASE SPACE DENSITIES 

In the previous sections, the phase space densities were 
evaluated at one particular time. For Hamiltonian dy- 
namics this is enough information to determine the aver- 
age dissipated work. In contrast, when the dynamics are 
stochastic the phase space densities alone do not contain 
enough information to determine the average dissipated 
work. In this section we address this issue by investigat- 
ing the implications of specifying the phase space density 
at two times. Let P[zi, ti; Zo,to) be the "two-time" prob- 
ability for the system to be at zq at time to and z\ at a 
later time t\. We will explicitly demonstrate that two- 
time distributions provide a better bound for the average 
dissipated work: (W diss ) > D(P f \\Pr) > D{p F \\p R ). 

For clarity we consider the overdamped limit as in 



-If 

0k v 



1 



-2fc(ti-to)/7 



(25b) 



The pairwise distribution is the product of Gaussian dis- 
tributions; as such, it is Gaussian as well. A simple, yet 
lengthy calculation shows that the means and variances 
of the pairwise distribution for the forward process are 



uto-^(l 
ut x - f (1 



e -*ii/7 



(26a) 



2 



1 / 1 e -fc(*i-*o)/7 

= ~0k V e-^* 1 -* ^ 1 

(26b) 

Repeating the above calculation for the reverse joint dis- 
tribution leads to similar results. 

The next step is to determine the relative entropy. 
We are interested in comparing the pairwise probabil- 
ity for the forward process Pf{x\, t%; Xq, to), with the 
time-reversed pairwise probability for the reverse process 
Pr(xo,t — to;xi,r — ti). If we allow to and t\ to be 
arbitrary, then from Eq. [TT] the relative entropy is 



D(P f \\Pr) = 



20j 2 u 2 /k 

I — e -2fc|ti-t |/7 



where 



x [al + a 2 



2 e ~ fc l t i-*o|/7 



a ai 



(27) 



= l-e- fe ^coshf-(l-t i )j, i = 0, 1. (28) 



Figure [3] compares the "two-time" relative entropy 
D(Pf\\Pr), "one-time" relative entropy D(pp\ \pr), and 
the average dissipated work. For D(P f \\Pr), we fix t = 
r/4. We see that 0{W diss ) > D(P F \\P R ) > D( Pf \\pb) for 
all times t\. Thus, by specifying an additional time in our 
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0.00 0.25 0.50 0.75 1.00 



is not true for stochastic systems. Calculating the dis- 
sipation from microscopic information in stochastic sys- 
tems requires knowledge of the entire evolution of the sys- 
tem, both for the forward and reverse processes. Partial 
knowledge gives only a lower bound on the dissipation. 
We have seen that the tightness of the lower bound is 
correlated with the amount of information known about 
the system's evolution. Specifically, when the system is 
loosely coupled to a reservoir, very little information is 
lost to the reservoir. As a result the relative entropy 
between the system's forward and reverse phase space 
densities reasonably approximates the dissipation in this 
limit. The lower bound can also be tightened using a 
multi-time relative entropy where the microscopic state 
of the system is specified many times along the system's 
trajectory. 

Acknowledgments 



FIG. 3: Comparison of D(Pf\\Pr) (dot-dashed line), 
D(pf\\pr) (solid line) and /3{Wdiss) (dashed line) in the over- 
damped limit. All system parameters have been set to one in 
their respective units. 



probability distributions we have improved the bound in 
Eq.ffl 

This result is an instance of the coarse-graining proce- 
dure outlined in Section [TT] (see the discussion below Eq. 

where the path space probability distribution T'f{i F ) 
has been projected down onto a two-time phase space 
density. Projecting onto a two-time phase space density 
eliminates less microscopic information than projecting 
onto a one-time phase space density; hence, the two-time 
relative entropy is greater than the one-time relative en- 
tropy. 

In a related independent analysis, Gomez-Marin et al. 
EH used the present model to investigate the behav- 
ior of the "TV-time" relative entropy Dm- By evaluating 
Dn at N equally spaced times, i.e. ti + i — t, = t/N, they 
were able to demonstrate that f3(Wdi SS ) — -Djv ~ l/N 2 
as N — ► oo. We note that our two-time relative entropy 
agrees with theirs: choosing t\ = r and to = in Eq. 
we recover Eq. 30 of Ref. [H and Eq. 22 (with n = 1) 
ofRef. 0. 



V. CONCLUSION 

The average energy dissipated by a thermodynamic 
process is related to how distinguishable the process is 
from its time-reversal, i.e. how well one can discern the 
"arrow of time" . This paper provides a pedagogical illus- 
tration of this idea using an exactly solvable model. For 
deterministic Hamiltonian systems, complete knowledge 
of the forward and reverse phase space density at any 
one time is enough to determine the dissipation. This 



We thank S. Vaikuntanathan and A. Ballard for a 
critical reading of this manuscript and the University of 
Maryland for their generous support. 



APPENDIX A: PHASE S PACE DENSITY FOR 

section Unci 

In this appendix, we solve Eq. [M] for the full phase 
space density of Section fill CI To begin, the coefficients 
of Eq. are linear in x and p. Consequently, as in the 
previous sections, an initial Gaussian distribution will 
remain Gaussian. Again we must determine the means, 
x and p, and covariance matrix a as functions of time. 
Combining Eq. [M] with the derivatives of Eqs. [5] we find 
the equations of motion for the means and variances of 
the forward process 

dxp = Pf^^ dpF = _ k ^- p _ ut ^ _ _7 ( Ala ) 
at m at m 



dt 



dt 



27 2 



da T 



7 



—al - 2ka xp + 

771 p 
1 



(Alb) 



a xv - kai H al. 



dt m to *" 

The reverse process phase space density is described by 
a similar set of equations. 

The solutions of Eqs. lAlal and their reverse process 
counterparts are 

x F (t) = Ae r+t + Be r - 1 +u(t-j/k), (A2a) 

p F (t) = mr + Ae r+t + mr~Be r - t + mu, (A2b) 

x R (t) = -Ae r + t -Be r - t + u(T-t + -//k),(A2c) 

p R (t) = -mr + Ae r+t - mr^Be r - t - mu, (A2d) 



where 



r± 



± -\/7 2 — 4777fc 

2m 



(A3) 



7 



and 



Using Eq. the average dissipated work is 



4mk 



r yr + /k + 1 
— 7r_/fc — 1 



(A4) 



are determined by the initial conditions (Eq. I13ap . The 
solutions to Eq. lAlbl with initial conditions given by Eq. 
Il3bl are 

4(t) = J, °l(t) = -k> cr xp (t)=0. (A5) 



6k 



(Wdiss) = iu 2 r + uk 
From Eq. 1111 the relative entropy is 



— (1 - e^ 4 ) + — (1 - e r 



(A6) 



D(pf\\pr) = \\pk [A (>+< + e r +^) + B (e r 



+ fim 



(A7) 
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